x{n) 
• o • 



I • I I — hH- 



-hH — I I I I — I I I I • T I — In 



-3-2-1 0 1 2 3 



• • • • 



y(t) 



y 3T -2T -IT 0 T 2T 3T f 




FIG. 1 




FIG. 2 



{ INTERPOLATE.dsp 

Real time Direct Form Filter, N taps, uses an efficient algorthm 

to interpolate by L for an increase of L times the input sample rate. A 

restriction on the number of taps that N/L be integer. 

INPUT: adc 
OUTPUT: dac 

! 

MODULE/RAM/ABS=0 interpolate; 



.CONST N=300; 

•CONST L=4: { interpolate by factor of L ] 

.CONST NoverL=75; 

.VAR/PM/RAM/CIRC coef[N]; 

.VAR/PM/RAM/CIRC data[NoverL]; 

.VAR/PM/RAM/ counter; 

.PORT adc; 

.PORT adc; 

•INIT coef: <coef.dat>; 



RTI; 
RTI; 
RTI; 

JUMP sample; 

initialize: IMASK=b#0000; 

ICNTL=b#01111; 
Sl=1; 

DM(counterO=SI; 

l4=^oef; 

L4=^coef; 

(listing contiunes on nest page) 



[interrupt 0 ] 

[interrupt 1 } 

[interrupt 2 ] 

[interrupt 3 at (L*input rate) \ 

[disable all interrupts^ 
[edge sensitive interrupts? 
[set interpolate counter to 1j 
[for first data sampiej 
[setup a circular buffer in PM| 



FIG. 



4 



M4=L; Smodifier for coef is L| 

M5=-1; jmodifier to shift coef back -]\ 

IO=^data; jsetup circular buffer in DMj 
LO=^data; 
MO-1 

IMASK=B#1000; Sencble interrupt 3l 
waitjnterrupt: JUMP waitJnterrupt;Sinfinte wait \oop\ 

\ Interpolate \ 

sample: M0DIFT(I4,M5); Jshifts coef pointer back by -1| 

AYO=DM(counter); 

AR=AY0-1; {decrement and update counter} 

Oy(counter)=AR; 

IF NE JUMP dojir; ftest ant input if L timesj 

\ input data sample, code executed at the sample rate ] 

dojnput: AYO=DM(adc); Jinput data sample^ 

DM(IIO,MO)=AYO; {update delay line wiht newestj 

M0DIFY(I4,M4); jshifts coef pointer up by Li 

DM(counter)=M4; {reset counter to l\ 

\ filter pass, occurs at L times the input sample rate \ 

do_fir: CNTR=NOVERL -1; \H/l since round on lost tap! 

MR=a MXO=DM(iaMO); MY0=PM(I4,M4); 

DO tapioop UNTIL CE; {N/L-1 taps of FIRj 
taploop: MR=MR+MXO*MXO(SS), MXO=DM(IO,MO), MY0=PM(I4,M4); 

IF MV SAT MR; {saturate result if overflowed! 

DM(dac)=MR1; {output sample! 

RTI; 

ENDMOD: 
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